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ABSTRACT 



Aims. In the interstellar clouds, molecular hydrogens are formed from atomic hydrogen on grain surfaces. An atomic hydrogen hops around 
till it finds another one with which it combines. This necessarily implies that the average recombination time, or equivalently, the effective 
7-H ■ grain surface area depends on the relative numbers of atomic hydrogen influx rate and the number of sites on the grain. Our aim is to discover 
^ 'this dependency. 

Methods. We perform a numerical simulation to study the recombination of hydrogen on grain surfaces in a variety of cloud conditions. We 
use a square lattice (with a periodic boundary condition) of various sizes on two types of grains, namely, amorphous carbon and olivine. 
Results. We find that the steady state results of our simulation match very well with those obtained from a simpler analytical consideration 
provided the 'effective' grain surface area is written as ~ 5", where, 5 is the actual physical grain area and a is a function of the flux of atomic 
\Q ' hydrogen which is determined from our simulation. We carry out the simulation for various astrophysically relevant accretion rates. For high 
accretion rates, small grains tend to become partly saturated with H and Hi and the subsequent accretion will be partly inhibited. For very low 
■ accretion rates, the number of sites to be swept before a molecular hydrogen can form is too large compared to the actual number of sites on 
the grain, implying that a is greater than unity. 
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1 . Introduction ditions and using rate and master equations presented by Biham 

et al. (2001). If one considers only H and H%, then, in a steady 



It is long known that the formation of H 2 in the interstellar 
cloud takes place on the grain surfaces and are then released 
into the gas phase. Considerable studies were made since then 
to understand the real physical processes which are taking place 



state, the number of hydrogen atoms («#) and molecules (nm) 

on a grain surface are obtained by equating the effective influx 

(or, accretion) rate of H with the rate at which they are used 

up to form H 2 and/or get desorbed from the grain surface into 
both theoretically (e.g., Hollenbach, Werner & Salpeter, 1971; ... " 

J ° ±- ' ' me „ as phase, i.e., 

Takahashi, Matsuda & Nagaoka, 1999; Biham et al. 2001) and 6 * 

experimentally (e.g., Pirronello et al. 1999). A quantitative esti- ff _ \Y H n H + 2a H n 2 H (la) 
mation based on physisorption of the formation rate of H 2 was 

carried out by Biham et al. (2001) and Green et al. (2001). It and 

was found that a significant production is possible in cooler WninHi = A^h^h- (lb) 
(~ 10 - 25K) clouds (Rae et al. 2003). Cazaux & Tielens 

(2002, 2004) use both physisorption and chemisorption, to Here ' is the effective flux (in units of number/s) of H on 

demonstrate that H 2 production is possible at high tempera- the g rain surface ' W n and W »2 (both in units of s 1 ) are the 

tures (~ 200 - 400K) also. Recently, Acharyya, Chakrabarti & desorption co-efficients of H and H 2 respectively, a H (in units 

Chakrabarti (2005) and Acharyya & Chakrabarti (2005) com- of s ~') is the effective rate of recombination of two H atoms 

puted the formation of H 2 from H assuming various cloud con- t0 form one H 2 molecule and 1 - y. is the fraction of H 2 that 

are formed on the grain which are spontaneously released to 

Send offprint requests to: Prof. Sandip K. Chakrabarti the gas phase. Solving these simple equations, we find in the 
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steady state, number of H and H 2 on a grain surface to be, 



^jw 2 H + 8a H <p H - W H 



nm = 



IM H n H 
W H 



(2a) 



(2b) 



We ignore H 2 formation by tunneling since the experiments 
indicate hopping to be the major process (see, Pironello et al. 
1999; Katz et al. 1999). In the literature, it is customary to de- 
fine the effective rate of recombination to be, 



a H =A H /S, 



(3a) 



(See, Acharyya, Chakrabarti, & Chakrabarti, 2005 and refer- 
ences therein), where A H is the hopping rate (inverse of the dif- 
fusion time Td — 1 /Ah) on a grain surface and S is the number 
of sites on the grain surface (Biham et al. 2001). The argument 
for reducing the diffusion rate by a factor of S is this: on an 
average, there are n — S number of sites in each direction of 
the grain. Since the hopping is random, it would take a square 
of this, namely, n 2 = S number of hopping to reach a distance 
located at n = S 1 ^ 2 sites away, where, on an average, another 
H is available. Thus, the effective recombination rate was cho- 
sen to be A H /S . However, when the flux is very high in a small 
grain, many sites would be full and an H need not hop a dis- 
tance of S to meet another one. Similarly, when the accretion 
rate is too low, one H may have to sweep the whole grain sev- 
eral times before it can meet another H to recombine. Thus, it 
is expected that the effective rate of recombination would be a 
function of the number of sites as well as the flux of H on the 
grain. If we assume that the effective number of sites could be 
written as S a °, where ao is a constant for a given flux and grain, 
then, in principle, we could define the effective recombination 
rate to be, 

a H =A H /S a \ (3b) 

and we expect that higher the accretion rate, lesser would be 
the value of ao- In the opposite limit, when the accretion rate 
is low, one would get ao > 1 since the effective site number 
is higher than S . In order to see this effect we clearly need to 
actually compute time dependences of H and H 2 by numerical 
simulation and derive ao when a 'steady state' is reached. This 
we do in this paper. 

In the next Section, we present the logic based on which the 
simulation is carried out. In §3, we present our results under 
various grain conditions and compare our results with the ex- 
pectation from the steady state (Eqs. 2ab). We show that these 
two results match only if ao is used in the steady state equa- 
tions. Finally, in Section §4 we draw our conclusions. 

2. Numerical Simulation 

The simulation is based on several physical processes which 
might be taking place on the grain surface. We assume each 
grain to be square in shape having S = n 2 number of sites 
(where, n is the number of sites in each direction). We use the 
periodic boundary condition. We consider only the mono-layer 



on the grain surface. Let A H = vexp(-Eo/kt,T), the rate at 
which H hops from a site to the next, where, v is the vibra- 
tional frequency, 



v = 



2sE, 



n 2 niH 



Here, s ~ 10 cm is the surface density of sites on a grain, 
m H is the mass (in gm) of the H atom and Ej is the binding en- 
ergy corresponding to physisorption (in meV). We choose the 
minimum time step to be 1 /Ah, the hopping time and advance 
the global time by this time step. Let Fh be the accretion rate 
of H which increases the number of H on a grain by sticking 
to it. Here, Fh - A < V > Nh s 4 , Nh is the number density 
(in cm 4 ) of H in the gas phase, < V > (in cm 4 ) is the average 
thermal velocity of the particle in the gas phase and A is the 
physical grain surface area (in cm 2 ); when 10 2 < N H £ 10 7 , 
and V ~ 10 4 , we have 3 x 10~ 8 < F H < 3 x 10~ 3 . We 
note that the rate at which the dropping of H takes place on 
a grain is a very small number. To handle this, we drop one 
H after every A H /F H number of hops. When F H is very low, 
we continued our simulation for a larger time interval till a 
steady state is reached. Similarly, where Fh is large enough 
so that Fh/Ah > 1, we drop several at each time step. In any 
case, for a given F H , the effective accretion rate is given by 
4>h = F H (1 - fgrh- fgrhi), where f grh and f grh2 , are the fractions 
of the grain which are occupied by H and H 2 respectively. After 
1/A H seconds, <fe/A H number of hydrogen atoms are dropped 
on the grain. The exact site at which one atom is dropped is 
obtained by a pair of random numbers (R x ,R y ; R x < 1 and 
R y < 1) obtained by a random number generator. This pair 
would place the incoming hydrogen at (i, j)Xh grid, where, i 
and j are the nearest integers obtained using the Int function: 
i = int(R x *n + 0.5) and j = int(R y * n + 0.5). Each atom starts 
hopping with equal probability in all four directions. This di- 
rection is decided by another random number. When, during 
the hopping process, one atom enters a site which is already oc- 
cupied by another atom, we assume that a molecule has formed 
and we increase the number of H 2 by unity and reduce the num- 
ber of H by two. However, when the atom enters a site occupied 
by an H 2 , another random number is generated to decide which 
one of the other nearest sites it is going to occupy. Thermal 
evaporation of H and H 2 from a grain surface are handled in the 
following way. Let Wh and Wm be the desorption co-efficients 
of hydrogen and hydrogen molecule given by vexp(-E\/kbT) 
and vexp(-E 2 /kbT), where E\ and E 2 are the activation barrier 
energy (in meV) for desorption of H and H 2 molecule. Thus 
one atom is released to the gas phase after every 1 / Wh seconds 
and one molecule is being released after every I/Whi seconds. 
We generate a random number R t for every H present on the 
grain and release (at each time step, i.e., 1/A H seconds) only 
those for which R, < Wh/A H - A similar procedure is followed 
for the evaporation of H 2 for which the criterion for evaporation 
was R, < Whi/Ah- There could also be spontaneous desorption 
and a factor of of n H2 is evaporated to the gas phase. Here 
too, a random number R s is generated for each newly formed 
H 2 present on the grain. Those which satisfy R s < 1 - /z are 
removed to the gas phase. 
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Fig. 1. (a-b): Variation of ao as a function of the accretion rate 
per site for an (a) amorphous carbon grain kept at 20K and (b) 
olivine grain kept at 10K. The solid, dot-dashed and dashed 
curves are for 10 4 , 9 x 10 4 and 10 6 sites respectively. As the 
accretion rate is reduced, ao starts to deviate from 1 as the ef- 
fective grain site number goes up due to sweeping of the grain 
more than once before recombination takes place. The devia- 
tion is highlighted using dotted curves by extrapolating at very 
low accretion rates. 

We continue our simulation for each type of grain and for 
each accretion rate and temperature for a duration more than 
10 6 s or more depending on the accretion rate Fh, till a quasi- 
steady state is reached. After an initial transient period, rin and 
nm, the number of H and H2 reach a quasi-steady state with 
some fluctuations. These numbers can be used to compute the 
recombination time. The values of the activation barrier ener- 
gies Eo, Ei and E 2 are taken from Katz et al. (1999) and are 
given by, £0 = 24.7 meV, E\ = 32.1 meV and E2 = 27.1 meV 
for olivine and Eq = 44 meV, E\ = 56.7 meV and £2 = 46.7 
meV for amorphous carbon grains. For olivine, /i = 0.33 and 
for amorphous carbon n = 0.413 was used. 

3. Recombination time scale as a function of grain 
and gas parameters 

Using the results of our simulation, we are now in a position to 
compute a quantity a(t) defined by, 



a(t) = log(- 



2A H n 2 H 



-)/log(S). 



(4) 



4>h - W H riH 

As t — » 00, a(t) — » ao. If there were no effect of the accre- 
tion rate and the finite size of the grain, then we would ob- 
tain a(t) = 1 identically at all the time. However, we find dif- 



Fig. 2. (a-b): Variation of (3q as a function of the accretion rate 
per site for an (a) amorphous carbon grain kept at 20K and (b) 
olivine grain kept at 10K. The solid, dot-dashed and dashed 
curves are for 10 4 , 9 x 10 4 and 10 6 sites respectively. As the ac- 
cretion rate is reduced, /3q becomes very high compared to unity 
for the same reason as ao- As the accretion rate is increased the 
grain becomes mostly occupied and an H has to hop, on an av- 
erage of S 1/2 times before an H2 is produced. Thus /?o becomes 
almost 1/2. The deviation is highlighted using dotted curves by 
extrapolating at very low accretion rates. 



ferently. In Figs, l(a-b) we plot ao for amorphous carbon and 
olivine grains at T = 20K and T = 10K respectively as a func- 
tion of the accretion rate per site a s . The solid, dot-dashed and 
dashed curves are for 10 4 , 9 x 10 4 and 10 6 sites respectively. 
For very very low rates, it becomes impossible to carry out the 
simulations in a computationally finite time. We thus extrapo- 
lated the curves (dotted) for very low a s values. We find that 
as the accretion rate per site goes down, ao deviates from unity 
significantly. Thus, clearly, for theoretical considerations, the 
number of H and H2 should be calculated using Eq. (3b) and 
not Eq. (3a). 

In this context, one important parameter /? may be defined 
as the 'catalytic capacity' of a grain which measures the effi- 
ciency of the formation of H2 on a grain surface for a given 
pair of H residing on it. Let 6Nh2 be the number of H2 formed 
in 8t time. Since two hydrogen atoms are required to create one 
H2, the average rate of creation of one H2 per pair of H atom 
would be given by, 



<A m >= 



1 6Nh 2 
2hh St 



(5a) 
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Fig. 3. : Comparison of the simulation results (dark circles) 
with those obtained from analytical considerations [Eq. 2(a- 
b)] (dashed curves) when suitable modification of the average 
recombination rate (Eq. 3b) is made. An olivine grain of 10 4 
sites at a temperature of 10K has been chosen in this compar- 
ison. Dotted curves are drawn using analytical results for ao 
extrapolated to very low accretion rates. 

We identify the inverse of this rate with the average formation 
rate given by, 

T f (t)=S m /A H . (5b) 

Thus, 

S fi (t) = A H I <A m >. (5c) 
This yields /3(f) as a function of time as, 

/3 = log(A H / <A m >)llog(S). (6) 

Here too, we can define /?o = yS(f — > oo). Our goal would be 
to see if fio thus obtained actually varies with grain parameters, 
accretion rate and temperatures. We anticipated in Acharyya, 
Chakrabarti and Chakrabarti (2005) that this exponent could 
be as low as 0.5 or less. Indeed, this is what we see as well. 
In Figs. 2(a-b) we show {3 Q as a function of the accretion rates 
per site for various grain sites and for (a) amorphous carbon at 
T = 20K and (b) olivine at T = 10K grains respectively. The 
solid, dotted and dashed curves are for 10 4 , 9 x 10 4 and 10 6 
sites respectively as before. We find that j3o can vary anywhere 
between ~ 1.5 to ~ 0.5. Most interestingly, as the accretion 
rate becomes high, the exponent becomes close to 0.5, as if the 
recombination is taking place on a one-dimensional grain! This 
is perhaps understandable, since as the grain gets filled up, the 
maximum number of hopping in any direction should not be 
more than n, i.e., the number of sites in any one direction of the 
grain. 

It is interesting to compare the results of our simulation 
with those obtained from the analytical considerations with 
and without our ao factor in the equation mentioned above. In 
Table 1, we present this comparison. We take an olivine grain 
of 10 4 sites at 10K and vary the accretion rates. In Column 1, 
we give the accretion rate per site of the grain. In Column 2 we 
present the coefficient ao which we derive from our simulation. 
In Columns 3-5, we present the number of H as obtained by our 
simulation and the modified equation (Eq. 3b) and the standard 
equations (Eq. 3a). Columns 6-8, we present similar results for 



H2- We find that our simulation matches more accurately with 
the analytical results provided Eq. (3b) is chosen. If standard 
equation (Eq. 3a) is used, the deviation is very significant. 

In Fig. 3, we compare the number of H and H2 on an olivine 
grain of 10 4 sites kept at 10 K as a function of the accretion rate 
per site. The dashed curves are drawn using Eqs. 2(a-b), and 
3b (using ao from the simulation) while the dark circles are 
purely from our simulation. Using extrapolated ao (from Fig. 
lb) we computed n H and nm analytically and see very good 
agreements with the simulation at lower rates also. 

4. Concluding remarks 

In this paper, we carried out a numerical simulation to show 
that the effective grain surface is not really the same as the 
physical grain surface. The simulations have been carried out 
for amorphous carbon and olivine grains kept at 20K and 10K 
respectively for a wide range of accretion rates. In the litera- 
ture, it is assumed that the recombination time is the grain site 
number divided by the diffusion rate (Eq. 3a). However, we 
find that this simplistic assumption is not valid especially when 
the accretion rate of H on the grains is very low. In this case, 
an atomic hydrogen may have to sweep a grain several times 
before meeting another atom to form a molecule thereby in- 
creasing the effective surface area. For very high rates, the grain 
surfaces could be partially filled with H and H2 and search for 
another H need not take S number of hops. We found that on 
an average, ao, the index which determines effective number 
of sites, could be ^ 1 for very low rates. We also find that ao 
depends on the nature of the grains, the temperature and the 
grain site numbers. Thus the recombination time is a complex 
function of these parameters. From the simulation results, we 
defined another index (3 which is a measure of the average cat- 
alytic ability to form an H2. This index is also found to be a 
strong function of the accretion rate. In fact, for very high rates 
the index may go down to 0.5 as anticipated in ACC05. A com- 
parison of the recombination efficiency as obtained from our 
procedure with that obtained from the standard consideration 
[Eq. l(a-b)] clearly shows a deviation easily attributable to the 
site/rate dependence of ao in a realistic case. The considera- 
tions based on which we came to our conclusion is so generic 
that we believe that similar dependence should follow even for 
other grain surface reactions. This will be discussed elsewhere. 
We believe that this important behaviour will have to be kept 
in mind in future studies of chemical evolution in interstellar 
clouds in presence of grains. 

Work of AD has been supported by a DST Project and that 
of KA is supported by an ISRO Project with Centre for Space 
Physics. 
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Table 1. Comparison of H and H 2 abundances in various methods 
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